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Abstract. I review the color dipole formulation of heavy quark production in 
the light of recent RHIC data. Since charm and bottom production directly 
probe the gluon density, these processes allow one to study shadowing and 
parton saturation at RHIC and LHC. The dipole approach provides a convenient 
framework to calculate these nuclear effects. I present numerical results for open 
charm and bottom production in proton-proton and proton-nucleus collisions and 
discuss transverse momentum broadening of heavy quarkonia. 



At high center of mass energies y/s, the cross section for any reaction a + N — > 
{b, c, . . .}X can be expressed as convolution of the light-cone (LC) wavefunction 
for the transition a — > {b, c, . . .} and the cross section for scattering the color 
neutral {anti— a, 6, c . . .}-system on the target nucleon N. In the case of heavy quark 
production this means that the Feynman graphs in Fig. D can be written in the form 



DEI, 



where a qq a is the cross section for scattering a color neutral quark-antiquark-gluon 
system on a nucleon 0, 



a(GN^{QQ}X)= / da / d 2 p \^ G ^ Q n(a, p)\ a qqG {a, p), (1) 



<JqqG{a, p) = - [(J q q{ap) + <T qq (ap)} - -<T qq {p). 



(2) 



Here a is the light-cone momentum fraction carried by the heavy quark Q, and 
a — 1 — a is the momentum fraction of the Q. The flavor independent dipole 
cross section <r qq {p) is a nonperturbative quantity and has to be determined from 
experimental data. It depends on the transverse separation p between quark and 
antiquark. Note that the dipole approach is formulated in the target rest frame with 




Figure 1. The three lowest order graphs contributing to heavy quark production 
in the dipole approach. These graphs correspond to the gluon-gluon fusion 
mechanism of heavy quark production in the parton model. 
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Figure 2. Results for the total heavy quark pair cross section as function 
of cm. energy. Varying free parameters in the dipole approach (solid lines) 
and in the NLO parton model (dashed lines) [3f gives rise to the uncertainties. 
While the curves have already been published in Ref. this figure now also 
includes new preliminary data points J^j from PHENIX and STAR (left panel, 
^fs = 200 GeV) and a preliminary HERA-B measurement of open bb production 
(the lower HERA-B point.) 

the longitudinal axis parallel to the projectile gluon. Higher order corrections cause a q q 
to depend on the target gluon momentum fraction X2 as well, but this is not explicitly 
written out here. The LC wavefunctions ^g-^qq can be- calculated perturbatively, 
explicit expressions can be found e.g. in Ref. [2]. 

In momentum space, Eq. Q can be written in fc^-factorized form. Hence the 
dipole formulation is closely related to the approach of Ref. [3J- However, the dipole 
approach takes into account only the finite transverse momentum of the target gluon. 
The relation of the dipole approach to the conventional parton model is explained in 
Ref. @j. 

In order to calculate the cross section for heavy quark pair production in proton- 
proton (pp) collisions, Eq. Q has to be weighted with the projectile gluon density. 
Numerical results are shown in Fig. [JJ Varying free parameters, such as the heavy 
quark mass mg , leads to the uncertainty represented by the space in between the solid 
curves (or in between the dashed curves in case of the next to leading order (NLO) 
parton model, respectively), see Ref. @] for a complete discussion. As already pointed 
out in Ref. 0|, HERA-B energy is apparently too low for the dipole approach. 

Partonic configurations with fixed transverse separations in impact parameter 
space are known as eigenstates of the interaction. Therefore, the dipole approach is 
especially suitable to describe multiple scattering effects, which become important at 
low xi , when the lifetime of the QQ-pair becomes large enough to allow for more than 
one interaction with the target nucleus. There are two different sources of nuclear 
suppression for heavy quark production [2J: 

First, the QQ pair in Fig. n can rescatter several times inside the target. Since 
the amplitude for interaction with a single nucleon is known and because dipoles 
are interaction eigenstates, all these rescatterings can be resummed by eikonalizing 
o-gg in Eq. J5J. Note that since the typical size of the pair is p 2 ~ I/tuq, double 
scattering behaves parametrically like Acr/rr cx Q 2 (x2)/rriQ. Though formally higher 
twist, nuclear suppression due to heavy quark rescattering is enhanced by a factor of 
the saturation scale in the target nucleus, Q 2 s {x2) oc A 1 / 3 and cannot be neglected for 
charm quarks. 
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In addition, there is the leading twist gluon shadowing One of the quarks 
(or the gluons) in Fig. 2] can radiate another gluon, which then propagates together 
with the QQ-pair through the nucleus. The magnitude of the rescattering correction 
is determined by the transverse distance pq this gluon can propagate from its parent 
quark. Unlike the size of the heavy quark pair, pa is limited only by nonperturbative 
QCD effects. In Ref. 7 , these nonperturbative effects were modeled for the case 
of DIS by an attractive interaction between the gluon and the quark it is radiated 
off. In order to explain the smallness of the Pomeron-proton cross section, a value 
Pq < 0.3 fm is needed, resulting in a rather weak gluon shadowing in DIS. However, 
the concept of a nonperturbative interaction suggests that pa may depend on the color 
state the QQ-pair is in after radiation of the gluon jHJ |5] . A pair in an octet state 
will interact more strongly with the gluon (small po), while the interaction would be 
much weaker between a color singlet pair and the gluon (large pq). That way, gluon 
shadowing would be much stronger for charmonium production in pA than in DIS. 
Here, I shall take a different point of view by arguing that pq is not determined by 
a nonperturbative interaction, but by properties of the QCD vacuum that prevent 
the gluon from propagating distances > 0.3 fm (independent of the color state of the 
heavy quark pair). Gluon shadowing is then process independent, and is calculated as 
explained in the appendix of Ref. [Hj ■ This viewpoint is supported by the observation 
of a rather weak suppression of the J/\I/ yield in dAu collisions [TJH- However, these 
data do not disprove the predictions of Ref. 8 , which are for \ mesons rather than 
for J/*. 

Gluon shadowing Rq{x2 1 o) at impact parameter b is included in the calculation 
through the modified eikonal approximation, 

<4(x 2j p)=2 J d 2 bS [ l-exp^cr qq (x 2: p)T(b)R G (x 2 ,b^jY (3) 

where T(b) is the nuclear thickness. Note that the leading twist gluon shadowing has 
the effect of making the target nucleus more dilute, thereby reducing the saturation 
scale Q 2 S and diminishing higher twist saturation effects. 

Results are shown in Fig.|Hl The eikonal approximation assumes that the lifetime 
of the heavy quark pair (the coherence length l c oc l/x 2 ) is much larger than the 
nuclear radius Ra- This is certainly the case at the LHC (y/s = 8.8 TeV). For 
Vcs = at RHIC, however, l c ~ Ra and a more sophisticated calculation using 
the Green's function technique developed in Ref. 1 1 1| would yield a slightly smaller 
suppression. Gluon shadowing sets in at smaller x 2 than quark shadowing, since the 
QQG-fluctuation lives much shorter than the QQ-pair alone (because of the larger 
invariant mass of the QQG-st&te \V2\.) This is taken into account in this calculation. 
In particular there is no gluon shadowing at y cS = at RHIC. The rescattering of the 
heavy quark pair is a sizable contribution to the nuclear suppression of open charm 
production, especially at y/s — 200 GeV, where gluon shadowing is weak. This is 
one reason why the observed suppression in J yields in dAu collisions at RHIC 
cannot be explained by a modification of the nuclear gluon distribution alone. The 
rescattering of 66-pairs at the LHC is expected to be only a small effect. 

The PHENIX collaboration also reported values for transverse momentum 
broadening of J/*f? mesons in dAu collisions ^Oj that are significantly larger than 
in fixed target experiments. An increase of 5(pj>) with energy is expected in the 
dipole approach jl.'il 114) but the relevant energy scale is ^Jx\s rather than the 
hadronic cm. energy s/s. For RHIC kinematics, ^Jx~[s is even smaller than in 800 GeV 
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Figure 3. Rapidity dependence of nuclear suppression for open heavy flavor 
production in pAu (yfs = 200 GeV) and in pPb (y/s = 8.8 TeV) collisions. 
Dashed curves show gluon shadowing only, while solid curves also include 
rescattering of the heavy guark pair. 



fixed target experiments. The predictions from Ref. |14|. assuming incoherent initial 
state rescattering of the projectile gluon, are 8{pj){x2 = 0.09) = 0.4 GeV 2 (exp: 
(1.77 ±0.35) GeV 2 ) and 5{p 2 T ){x 2 = 0.003) = 0.8 GeV 2 (exp: (1.29 ± 0.35) GeV 2 ). It 
is unclear what causes the large discrepancy at x 2 — 0.09. At smaller x 2 , there could 
be an additional contribution to 8(pj) from the color filtering mechanism of Ref. 9 . 
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